Dressed Polyakov loop and flavor dependent phase transitions 
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The chiral condensate and dressed Polyakov loop at finite temperature and density have been 
investigated in the framework of Nj = 2 + 1 Nambu-Jona-Lasinio (NJL) model with two degenerate 
u, d quarks and one strange quark. In the case of explicit chiral symmetry breaking with physical 
quark masses, it is found that the phase transitions for light u, d quarks and s quark are sequentially 
happened, and the separation between the transition lines for different flavors becomes wider and 
wider with the increase of baryon density. For each flavor, the pseudo-critical temperatures for 
chiral condensate and dressed Polyakov loop differ in a narrow transition range in the lower baryon 
density region, and the two transitions coincide in the higher baryon density region. 

PACS numbers: 12.38.Aw, 12.38.Mh, ll.30.Rd 



I. INTRODUCTION 

QCD vacuum is characterized by spontaneous chiral 
symmetry breaking and color confinement. The dynam- 
ical chiral symmetry breaking is due to a non-vanishing 
quark anti-quark condensate, {qq) ~ (250MeV)'^ in the 
vacuum, which induces the presence of the light Nambu- 
Goldstone particles, the pions and kaons in the hadron 
spectrum. The confinement represents that only colorless 
states are observed in the spectrum, which is commonly 
described by the linearly rising potential between two 
heavy quarks at large distances, Vqq{R) = (JsR, where 
tJs — (425MeV)^ is the string tension. 

It is expected that chiral symmetry can be restored 
and color degrees of freedom can be freed at high tem- 
perature and/or density. The interplay between chiral 
and deconfinement phase transitions at finite tempera- 
ture and density are of continuous interests for studying 
the QCD phase diagram [H-Q. The chiral restoration 
is characterized by the restoration of chiral symmetry 
and the deconfinement phase transition is characterized 
by the breaking of center symmetry, which are only well 
defined in two extreme quark mass limits, respectively. 
In the chiral limit when the current quark mass is zero 
m = 0, the chiral condensate {qq) is the order parameter 
for the chiral phase transition. When the current quark 
mass goes to infinity m ^ oo, QCD becomes pure gauge 
SU{i) theory, which is center symmetric in the vacuum, 
and the usually used order parameter is the Polyakov 
loop expectation value (P) [Ij], which is related to the 
heavy quark free energy. At zero density and chiral limit, 
lattice QCD results show that the chiral and deconfine- 
ment phase transitions occur at the same critical tem- 
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perature, e.g, see Ref. [15l - ll9f . and also review papers 
[23,[2l|. This result is highly nontrivial because these two 
distinct phase transitions involve different mechanisms at 
different energy scales. It has been largely believed for a 
long time that chiral symmetry restoration always coin- 
cides with deconfinement phase transition in the whole 
(T, ^) plane. 

However, for the case of finite physical quark mass, 
neither the chiral condensate nor the Polyakov loop is 
a good order parameter. For heavy quark, there is no 
dynamical chiral symmetry breaking (e.g, see thus 
no chiral restoration. On the other hand, the linear po- 
tential description for confinement property is not suit- 
able for light quark system. In recent years, several lat- 
tice groups have made much effort on investigating the 
chiral and deconfinement phase transition temperatures 
with almost physical quark masses, e.g, RBC-Bielefeld 
group 23], which later merged with part of the MILC 
group [22| and formed the hotQCD group [H, , and 
Wuppertal-Budapest group |29[ 4 3 3l . The result from the 
RBC-Bielefeld group in 2006 O found that the two 
pseudo-critical temperatures for Nj = 2 + 1 coincide 
at Tc = 192(7)(4)MeV. The Wuppetal-Budapest group 
found that for the case of iV/ = 2+1, there are three pseu- 
docritical temperatures, the transition temperature for 
chiral restoration of u, d quarks Tc^"'^^ = 151 (3) (3) MeV, 
the transition temperature for s quark number suscepti- 
bility = 175(2)(4)MeV and the deconfinement transi- 
tion temperature — 176(3) (4)MeV from the Polyakov 
loop. Recently, it is shown in pTi. [28j. by using an im- 
proved HISQ action, hotQCD collaboration results are 
close to the Wuppetal-Budapest collaboration results. 

The relation between the chiral and deconfinement 
phase transitions has also attracted more interest re- 
cently in studying the phase diagram at high baryon den- 
sity region [1(| . It is conjectured in Ref. [Il| that in large 
Nc limit, a confined but chiral symmetric phase, which is 
called quarkyonic phase can exist in the high baryon den- 
sity region. The quarkyonic phase or chiral density wave 
state is due to the quark-hole pairing near the Fermi sur- 
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face. Nevertheless, it attracts a lot of interests to study 
whether such a chiral symmetric but confined phase can 
survive in real QCD phase diagram, and how it competes 
with nuclear matter and the color superconducting phase 
[13,]. 

In the framework of QCD effective models, there is still 
no dynamical model which can describe the chiral sym- 
metry breaking and confinement simultaneously. The 
main difficulty of effective QCD model to include con- 
finement mechanism lies in that it is difficult to calcu- 
late the Polyakov loop analytically. Currently, the pop- 
ular models used to investigate the chiral and deconfine- 
ment phase transitions are the Polyakov Nambu-Jona- 
Lasinio model (PNJL] 34-41] and Polyakov linear sigma 
model (PLSM) [i^, [43j]. However, the shortcoming of 
these models is that the temperature dependence of the 
Polyakov-loop potential is put in by hand from lattice 
result, which cannot be self-consistently extended to fi- 
nite barron density. Recently, efforts have been made 
in Ref. [ij, to derive a low-energy effective theory for 
confinement-deconfinement and chiral-symmetry break- 
ing/restoration. 

Recent investigation revealed that quark propagator, 
heat kernels can also act as an order parameter as they 
transform non trivially under the center transformation 
related to deconfinement transition (46l - l48| . The excit- 
ing result is the behavior of spectral sum of the Dirac 
operator under center transformation [U l49l - l5l| . A new 
order parameter, called dressed Polyakov loop has been 
defined which can be represented as a spectral sum of 
the Dirac operator [5l|. It has been found the infrared 
part of the spectrum particularly plays a leading role 
in confinement (47l |. This result is encouraging since it 
gives a hope to relate the chiral phase transition with the 
confinement-deconfinement phase transition. The order 
parameter for chiral phase transition is related to the 
spectral density of the Dirac operator through Banks- 
Casher relation Therefore, both the dressed Polyakov 
loop and the chiral condensate are related to the spec- 
tral sum of the Dirac operator. Behavior of the dressed 
Polyakov loop is mainly studied in the framework of Lat- 
tice gauge theory |52h54| . Apart from that, studies based 
on Dyson-Schwinger equations (sB - fs?! and PNJL model 
[H,!!^ have been carried out. In those studies the role of 
dressed Polyakov loop as an order parameter is discussed 
at zero chemical potential. The dressed Polyakov loop 
at finite temperature and density has been investigated 
in the two-flavor NJL model in Ref. (60j . In this paper, 
we show the phase diagram in the framework of three- 
fiavor NJL model by using the dressed Polyakov loop as 
an equivalent order parameter. 

This paper is organized as follows: We introduce the 
dressed Polyakov loop as an equivalent order parameter 
of confinement deconfinement phase transition and the 
NJL model in Sec. |lll Then in Sec lIIIl we show the re- 
sults of three-flavor QCD phase diagram in T — fj, plane 
in the chiral limit and in the case of explicit chiral sym- 
metry breaking, respectively. At the end, we give the 



conclusion and discussion. 



II. DRESSED POLYAKOV LOOP AND THE 
THREE-FLAVOR NJL MODEL 

We firstly introduce the dressed Polyakov loop. To do 
this we have to consider a U{1) valued boundary con- 
dition for the fermionic fields in the temporal direction 
instead of the canonical choice of anti-periodic boundary 
condition. 



(1) 



where < < 27r is the phase angle and /3 is the inverse 
temperature. 

Dual quark condensate E„ is then defined by the 
Fourier transform (w.r.t the phase 0) of the general 
boundary condition dependent quark condensate |5ll453l | , 
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where n is the winding number. 

Particular case of n = 1 is called the dressed Polyakov 
loop which transforms in the same way as the conven- 
tional thin Polyakov loop under the center symmetry and 
hence is an order parameter for the deconfinement tran- 
sition ]^51-53]. It reduces to the thin Polyakov loop and 
to the dual of the conventional chiral condensate in infi- 
nite and zero quark mass limits respectively, i.e., in the 
chiral limit m — >■ we get the dual of the conventional 
chiral condensate and in the m — > cxd limit we have thin 
Polyakov loop [Ull. 

The Lagrangian of three-flavor NJL model [6l| is given 

as 

a 

- x{Det/[v;(H-75)^]+Det/[7^(l-75)V]}. (3) 

Where ip ~ {u, d, s)'^ denotes the transpose of the quark 
field, and m = Diag(mij, m^, rris) is the corresponding 
mass matrix in the flavor space, with a = 1, • • • , N'j — 1 
are the eight Cell-Mann matrices, and Det/ means de- 
terminant in flavor space. The last term is the standard 
form of the 't Hooft interaction, which is invariant under 
SU{3)l X SU{3)r X U{1)b symmetry, but breaks down 
the J7a(1) symmetry. 

The (/) dependent thermodynamic potential in the 
mean field level is given as following: 

f f 
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with 
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-2N, 



(27r)3 
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Where the sum is in the flavor space, Epj ~ + ^ 
- ^ _ _L r,. I „vj, _ ^i^i^ i-jjg constituent 



and i;^^ = Epj ± + i( 
quark mass 



M0,j = nil - 4G's(cr)0,i + 2K{a)^j{a)cf,,k, 
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where (i,j,k) is the quark flavor indices {u,d,s), and 
= {i)fij}f)cj,. We will only consider isospin symmet- 
ric quark matter and define a uniform chemical potential 
/i for u, d and s. 

It is known that the NJL model lacks of confinement 
and the gluon dynamics is encoded in a static coupling 
constant for four point contact interaction. However, as- 
suming that we can read the information of confinement 
from the dual chiral condensate, it would be interesting 
to see the behavior of the dressed Polyakov loop in a sce- 
nario without any explicit mechanism for confinement. 

The thermodynamic potential contains imaginary 
part. We take only the real part of the potential and 
the imaginary phase factor is not considered in this work. 
The mean field ((t}^ is obtained by minimizing the po- 
tential for each value of (/) € [0, 27r) for fixed T and fi. 
The conventional chiral condensate is {ct)tx — (V'V')ir- For 
brevity from here onwards we will represent the conven- 
tional chiral condensate as < u > . The dressed Polyakov 
loop El is obtained by integrating over the angle. 



III. PHASE DIAGRAM FOR THREE FLAVORS 

We investigate phase transitions for two cases, i.e., in 
the chiral limit and in the case of explicit chiral symmetry 
breaking with physical quark mass, and the correspond- 
ing parameters are taken from Ref.jG^] and [63l |64|: 





m,[MeV] 


m,[MeV] 






chiral-limit 








1.926 


12.36 


physical mass 


5.5 


140.7 


1.835 


12.36 



TABLE L Two sets of parameters in 3-flavor NJL model: the 
current quark mass niq for up and down quark and ms for 
strange quark, coupling constants G and K, with a spatial 
momentum cutoff A = 602.3 MeV. 
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FIG. 1: The conventional chiral condensate — {(j)u,d,a and the 
dressed Polyakov loop J]"-'*''' of u, d, s quarks as functions of 
temperature for different values of the chemical potentials. 
Here, — (cr) and Ei both are measured in [GeV^]. 




FIG. 2: Three-flavor phase diagram in the T — fi plane for the 
case of chiral limit. The solid line is the critical line for Ei, 
and the dashed line is the critical line for conventional chiral 
phase transition. 



A. Phase diagram in the chiral limit 

We firstly consider the case of chiral limit, i.e. = 
iTT-d = TOs = 0. In Fig. [U we show the behavior of 
the conventional chiral condensate —{a) and the corre- 
sponding dressed Polyakov loop Ei for u, d and s quarks 
at different chemical potentials as functions of tempera- 
ture. For both order parameters, it is observed there are 
three temperature regions for — (cr) and Si. For — (cr), 
at smaller temperatures it remains constant at a value 
corresponding to the value of the conventional chiral con- 
densate in the vacuum, then it drops to zero at the criti- 
cal temperature Tc, and eventually keeps zero above the 



critical temperature. The critical temperature decreases 
with the increase of the chemical potential. It is noticed 
that, in order to guide eyes, we have connected the two 
end-points of the order parameter at the jump. 

On the other hand the behavior for the dressed 
Polyakov loop is just the opposite. It remains zero for 
small temperatures and then jumps at the critical tem- 
perature, and finally saturates to a high value which 
varies very slowly with temperatures. The almost zero 
value of El for small temperatures is due to the fact that 
the U{1) boundary condition dependent general quark 
condensate nearly does not vary with the angle (j) for 
small temperatures (see Eq. [5]). 
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It is seen that the phase transitions for chiral restora- 
tion and dressed Polyakov loop are of 1st order in the 
whole T — fi plane. For two- flavor case, it was found 
these two phase transitions are of second order. The Nf 
dependent result is in agreement with the results given by 
Pisarski and Wilczek in Ref. [HI]. The first order phase 
transition in three-flavor case is due to the 't Hooft in- 
teraction in Eq. which contributes a cubic term in 
the thermodynamical potential in Eq. (j4]). 

Fig. [2] shows the phase diagram of three-flavor in the 
chiral limit. We find almost exact matching for the tran- 
sition temperatures calculated from these two quantities 
in the whole T — ^ plane. 



B. Phase diagram with physical quark mass 
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For the case of finite quark mass m„ = rud = 5.5MeV 
and rus = 140.7MeV, we have chosen the model parame- 
ters of GsA2 = 1.835, Kh^ = 12.36 with A = 602.3MeV 
as in Ref. [H to fit = 135.0MeV, 92.4MeV, 
rriK = 497.7MeV and m^, = 957.8MeV. 

In Fig. [3] and [H we show the behavior of the conven- 
tional chiral condensate —{a) and the dressed Polyakov 
loop Si at different chemical potentials as functions of 
temperature for u, d and s quarks, respectively. 
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FIG. 3: The conventional chiral condensate ~{a-)u,d and the 
dressed Polyakov loop E"''' of u, d quarks as functions of tem- 
perature for different values of the chemical potentials. Here, 
— (cr) and Ei both are measured in [GeV^]. 

For both cases, it is observed that there are three tem- 
perature regions for — (cr) and Ei. For — (cr), at smaller 
temperatures it remains constant at a value correspond- 
ing to the value of the conventional chiral condensate in 
the vacuum, then it rapidly decreases in a small win- 
dow of temperature and eventually almost saturates to 
a lower value. The decreasing occurs at different tem- 
peratures for different values of the chemical potentials. 
On the other hand the behavior for the dressed Polyakov 



FIG. 4: The conventional chiral condensate —{cr)a and the 
dressed Polyakov loop Ef of s quark as functions of temper- 
ature for different values of the chemical potentials. Here, 
— {a) and Ei both are measured in [GeV^]. 



loop is just the opposite. It remains almost zero for small 
temperatures and then rises rapidly, finally saturates to 
a high value which varies very slowly with temperatures. 
The almost zero value of Si for small temperatures is due 
to the fact that the U{1) boundary condition dependent 
general quark condensate nearly does not vary with the 
angle (p for small temperatures (see Eq. [5]). 

The critical temperature for a real phase transition 
or the pseudo-critical temperature for a crossover is ex- 
tracted from the susceptibility of the order parameter or 
the temperature derivative of the order parameter. For 
example, for chiral phase transition of strange quark, the 
(pseudo)critical temperature is extracted from the tem- 
perature derivative d{— < as >)/{dT). This quantity 
describes how fast the order parameter changes with tem- 
perature. Normally the critical temperature corresponds 
to the fastest change of the order parameter, and the tem- 
perature derivative of the order parameter shows a peak 
at the critical point. However, there are some subtleties 
to determine the pseudo-critical temperature for the chi- 
ral restoration of the strange quark. We show how we de- 
termine the pseudo-critical temperature of the crossover 
by using FiglSJ which is the temperature derivative of the 
chiral condensate of the strange quark corresponding to 
Figgl 

For /i = 0, from Fig [S] one can observe that the temper- 
ature derivative of the strange quark condensate shows a 
peak at r = 196MeV, correspondingly, from Fig|4l one 
can see that the strange quark condensate changes fast 
at r = 196MeV, which is the critical temperature for 
chiral phase transition of the u,d quarks at zero chem- 
ical potential. However, the value of the strange quark 
condensate at T";^ = 196MeV is still around its vac- 
uum value, one cannot locate the pseudo-critical tem- 
perature of the strange quark at T"^ = 196MeV even 
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though there is a peak for the temperature derivative of 
the strange quark condensate. The reasonable explana- 
tion of the fast change of the strange quark condensate 
at T^'-^ = 196MeV is that the strange quark feels the 
chiral phase transition of u, d quarks due to the flavor 
mixing effect. For /i = 0, from FiglS] one can also ob- 
serve a bump region of the temperature derivative of the 
strange quark condensate around T = 250MeV, however, 
there is no obvious peak shown up. Therefore, we cannot 
extract an explicit pseudo-critical temperature from the 
chiral phase transition of the strange quark. Correspond- 
ingly, we find that the strange quark condensate at /i = 
changes smoothly with temperature. 
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FIG. 5: The derivative of strange chiral condensate 
d{—{a)a)/dT as functions of T for different values of /x. 



in order to determine the pseudo-critical temperature of 
the crossover. This method is also used to determine the 
dressed Polyakov loop of the strange quark. The critical 
and pseudo-critical temperatures extracted from the tem- 
perature derivative of the order parameters are shown in 
Fig. [51 It is found that the the chiral and deconfinement 
phase transitions are flavor dependent. 

At low baryon chemical potential region when /i < 
270MeV, for light flavors, i.e. for u, d quarks, we ob- 
serve from Fig. [3] that the conventional chiral condensate 
and the dressed Polyakov loop change rapidly with the 
increase of temperature. From the temperature deriva- 
tive of the order parameters of the chiral condensate 
and dressed Polyakov loop, we can obtain two separate 
pseudo-critical temperatures and for fixed /i, and 
we find is always smaller than . 

However, in the chemical potential region when /i < 
270MeV, for s quark, from Fig. 2] we can see that the 
conventional chiral condensate and dressed Polyakov loop 
change smoothly with the increase of temperature. From 
the temperature derivative of the order parameters, one 
cannot extract the values of the pseudo-critical temper- 
atures as already discussed. Therefore, in Fig. [HI of the 
three-flavor phase diagram, we can read that in the re- 
gion around < < 270MeV, the phase transitions for 
u, d are crossover, and different order parameters have 
different pseudo-critical temperatures. The s flavor expe- 
riences a rapid crossover, and no pseudo-critical temper- 
atures can be extracted from the order parameters. From 
the lattice results in Ref.|31j] at zero chemical potential, 
there is also no pseudo-critical temperature for the order 
parameter of strange quark's chiral condensate. 



The temperature derivative of the strange quark con- 
densate at /J, = 200MeV is similar to the case at /i = 0. 
The only difference is that the peak moves to a lower 
temperature. The small jump at the large strange chi- 
ral condensate region is induced by the u, d quark chiral 
phase transition. It cannot be regarded as the phase tran- 
sition for strange quark even though it corresponds to a 
peak of the strange chiral susceptibility, because the or- 
der parameter does not change so much comparing with 
its vacuum value. It should still be regarded as in the chi- 
ral symmetry breaking phase. At /i = 320MeV, it is seen 
from Figinithat the left peak develops to a sharp peak at 
T"^'', and an obvious peak shows up in the right bump 
region. Therefore, one can extract the pseudo-critical 
temperature for the chiral phase transition of the strange 
quark. For higher chemical potential, e.g, /j, = 460MeV 
or ji — 490MeV, because u, d quarks are already in chi- 
ral symmetric phase, there is only one peak shows up for 
the temperature derivative of the strange quark conden- 
sate in Fig. [5l and the location of the peak gives the 
pseudo-critical temperature of the phase transition. 

As we have discussed in detail above, one has to com- 
bine the information from the order parameter itself as 
well as the temperature derivative of the order parameter 
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FIG. 6: Three-flavor phase diagram in the T — fi plane for the 
case of ruu = ma = 5MeV and rUs = 140.7MeV. The dash- 
dotted lines are the critical line for Ei, and the dashed lines 
are the critical line for conventional chiral phase transition in 
the region of crossover. The solid lines indicates the 1st order 
phase transitions, and the solid circle indicates the critical 
end points for chiral phase transitions of m, d quarks. 

At higher baryon chemical potential region, it is ob- 
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served from Figl3]that the conventional chiral condensate 
and the dressed Polyakov loop change sharply with the 
increase of temperature. From the temperature deriva- 
tive of the order parameters, we find that the phase tran- 
sitions are of first order, and the critical temperatures for 
chiral and dressed Polyakov loop coincide with each other 
around CEP. 

For s quark, from Fig. 3] we can see that when the 
chemical potential becomes higher and higher, the con- 
ventional chiral condensate and dressed Polyakov loop 
change more rapidly with the increase of temperature. 
The temperature derivative of the order parameters give 
separate values of the pseudo-critical temperatures in the 
region 270 < /i < 450MeV, and the two pseudo-critical 
temperatures merge in the region of > 450MeV. 

From Fig. |5] of the three-flavor phase diagram, we 
can read the critical end point for u, d flavors lies at 
{T^ip,lJicip) = (68.4MeV,317.8MeV), which is differ- 
ent from the results in Ref.^6fl] for pure two-flavor NJL 
model. The difference comes from: 1) different model 
parameters have been used, 2) the coupling of s quark 
io u,d quark contributes one extra term in the thermo- 
dynamical potential comparing with the pure two-flavor 
case. The location of CEP in this work is in good agree- 
ment with that in Ref. [66| . 

In FiglT] and FiglHl we show the details of locating 
the CEP. In the first order phase region, there are two 
branches of number densities, i.e, for fixed chemical po- 
tential, the number density = has a jump at 
the transition temperature. The two branches of num- 
ber densities merge at the CEP. This feature is shown in 
FiglT] At the CEP, the phase transition is of second or- 
der and this is indicated by the divergent behavior of the 
number susceptibility. We show the number susceptibil- 
ity Xq — functions of the temperature in Fig IS] 
It is clearly seen that Xq develops a sharp peak at CEP. 
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FIG. 7: The quark number density as functions of the 
temperature for different chemical potentials, and Uq is in 
unit of [GeV^]. 
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FIG. 8: The number susceptibility Xq/T'^ ^ functions of the 
temperature for different chemical potentials. 



IV. CONCLUSION AND DISCUSSION 




FIG. 9: Conjectured 3D QCD phase diagram at finite tem- 
perature T, quark chemical potential jiq and isospin chemical 
potential ni. 

We investigate the chiral condensate and the dressed 
Polyakov loop or dual chiral condensate at finite tempera- 
ture and density in the three-flavor Nambu-Jona-Lasinio 
model. It is found that in the chiral limit, the phase tran- 
sitions are of 1st order and the critical temperature for 
chiral phase transition coincides with that of the dressed 
Polyakov loop. In the case of explicit chiral symmetry 
breaking, it is found that the phase transitions are fla- 
vor dependent, and there is a phase transition range for 
each flavor. The transition range of s quark is located at 
higher temperature and higher baryon density than that 
of u, d quarks. At low baryon density region, it is found 
that the transition range of u, d quarks are not separated 
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too much from that of the s quark, however, the separa- 
tion of the transition ranges for u, d quarks and s quark 
become wider and wider with the increase of the chemical 
potential. 

For light u, d quarks, the pseudo-critical temperature 
for chiral transition is smaller than that of the dressed 
Polyakov loop in the low baryon density region where 
the transition is a crossover, and these two phase tran- 
sitions coincide in the 1st order phase transition region 
at high baryon density. For s quark, both transitions are 
of smooth crossover at low baryon density, and becomes 
rapid crossover at moderate baryon density region where 
the pseudo-critical temperatures for the chiral conden- 
sate and the dressed Polyakov loop are separated, then 
at enough high baryon density, these two transitions co- 
incide with each other. 

Our results are based on the NJL model, where the 
gluon dynamics is encoded in a static coupling constant 
for four point contact interaction, a quantitative com- 
parison will not match with lattice results. However, we 
believe the scenario of the sequential phase transitions is 
physically correct. 

Till now, there are six quark flavors observed in exper- 
iment. These six flavors cover a very wide energy scale, 
from several MeV to several hundred GeV. Only light 
quarks experience dynamically chiral symmetry breaking 
in the vacuum, and chiral phase transition in high tem- 
perature and density. However, there is no good order 
parameters to describe the deconfinement phase transi- 
tion of light quarks. The conventional Polyakov loop is 
a good order parameter for confinement deconfinement 
phase transition in the limit of infinity heavy quark mass, 
and has the interpretation of the free energy of an infin- 
ity heavy quark. In analogy to that we can regard the 
dressed Polyakov loop as an order parameter for confine- 
ment deconfinement phase transition for a quark with 
mass TO, and interpret the dressed Polyakov loop as the 
free energy of a quark with any mass to .54j . There- 
fore, in principle, each flavor can have different critical 
temperatures for deconfinement phase transition. Lattice 
results already reflect such properties at zero chemical 
potential, e.g. the pseudo-critical temperatures for order 
parameters of u, d quarks, s quark and the Polyakov loop 
are different, and the the pseudo-critical temperature is 
higher for heavier quark mass. 

It is natural to understand that the separation of 
the phase transition range for different flavors becomes 
wider and wider with the increase of the chemical po- 
tential. Lattice result at zero chemical potential gives 
that the pseudo-critical temperature for u, d quarks is 
around 155MeV, and for s quark is around 175MeV. 
The difference is around 20MeV. However, at zero tem- 
perature, the u, d quarks restores chiral symmetry at 
the chemical potential around their vacuum constituent 
masses, i.e. /x"''' ^ Mu^d ^ 330MeV, and the s quark re- 
stores chiral symmetry at the chemical potential around 
^i^r^ M,r^ 550MeV. The difference is around 200MeV. 



Based on above analysis, in Fig. |9l we show our con- 
jectured 3 dimension (3D) QCD phase diagram for flnite 
temperature T, quark chemical potential jiq and isospin 
chemical potential /i/. 

In the plane of (/i, T), each flavor has its own transition 
range. The transition range is wider in the low baryon 
density, and becomes narrower and narrower with the in- 
crease of the chemical potential, and eventually merge at 
higher chemical potential. By using the lattice results 
at zero density, we identify the phase transition range 
around 155MeV for u, d quark, 175MeV for s quark, and 
190MeV for heavy flavor. The upper solid line is for the 
Polyakov loop, which does not change so much with the 
increase of baryon density. This result agrees with that in 
any Polyakov loop NJL model and Polyakov loop linear 
sigma model. Due to the flavor dependent phase tran- 
sitions, we naturally expect the color superconducting 
phase for two-flavor quark system and three-flavor quark 
system in different baryon density regions [Tsj . Due to 
the flnite mass of strange quark, the three-flavor color 
superconducting phase can be in the color flavor lock- 
ing (CFL) phase |67i|. CFL-kaon condensate phase (CFL- 
K)@, or uSC/dSC phase [6|. 

When isospin asymmetry is considered, the phase di- 
agram becomes much more complicated. At low baryon 
density region, there will be pion superfluidity and kaon 
superfluidity phases [t^. In the color superconducting 
phase, because isospin asymmetry induces mismatch be- 
tween the pairing quarks, there will appear unstable gap- 
less excitations [7l[ [t^ when charge neutrality condition 
is considered. It has been vastly discussed in many liter- 
atures about the true ground state of the charge neutral 
two-flavor and three-flavor cold quark matter, e.g, the 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) sate or other 
crystalline structure [zl] , the gluon condensate state [tI] , 
the current generation state |75ll. and so on. The de- 
tailed analysis given in Ref. |76| show that in the gapless 
color superconducting phase, both the phase part and 
magnitude part of the order parameter will develop in- 
stabilities. The phase part develops into the chromomag- 
netic instability, which induces the plane- wave state; The 
magnitude part develops the Sarma instability and Higgs 
instability, the Sarma instability can be competed with 
charge neutrality condition. If the Higgs instability can- 
not be cured by the electric or color Coulomb interaction, 
it will induce the inhomogeneous state. 
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